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Abstract 

An argument by Banks, Susskind and Peskin (BSP), according to which 
violation of unitarity would violate either locality or energy-momentum con¬ 
servation, is widely believed to be a strong argument against non-unitarity of 
Hawking radiation. We find that the whole BSP argument rests on the cru¬ 
cial assumption that the Hamiltonian is not highly degenerate, and point out 
that this assumption is not satisfied for systems with many degrees of free¬ 
dom. Using Lindblad equation, we show that high degeneracy of the Hamil¬ 
tonian allows local non-unitary evolution without violating energy-momentum 
conservation. Moreover, since energy-momentum is the source of gravity, we 
argue that energy-momentum is necessarily conserved for a large class of non- 
unitary systems with gravity. Finally, we explicitly calculate the Lindblad op¬ 
erators for non-unitary Hawking radiation and show that they conserve energy- 
momentum. 


PACS Numbers: 04.70.Dy, 03.65.Yz 


1 Introduction 

According to the semi-classical theory, Hawking radiation from black holes [1] is a 
non-unitary process [2]. Nevertheless, there are many different attempts to restore 
unitarity of Hawking radiation (see e.g. [3, 4, 5, 6, 7, 8, 9, 10, 11, 12] for reviews) and 
the belief that the true black-hole evolution should be unitary is highly prevailing. 
Yet, neither of the attempts to restore unitarity is without intrinsic problems, and 
many of the attempts contradict each other. In a confrontation with these problems, 
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it seems wise not yet to completely dismiss the possibility that Hawking radiation 
might really be a non-unitary process. 

One of the most frequent arguments for unitarity is authoritative reference to 
the work by Banks, Susskind and Peskin (BSP) [13], who argued that violation of 
unitarity would violate either locality or energy-momentum conservation. Some pos¬ 
sible loopholes in BSP arguments have been suggested in the literature [14, 15, 16, 
17, 18, 19, 20], but neither of these loopholes seems to be taken seriously in a wider 
community. 

In this paper we find and exploit a new loophole in the BSP argument. In Sec. 2 we 
critically re-examine the arguments used in [13] and find that one crucial assumption 
in their argument is not satisfied. This is the assumption that the Hamiltonian is not 
highly degenerate, which we show in Sec. 3 to be a wrong assumption. In Sec. 4 we 
explain in detail why high degeneracy of the Hamiltonian allows non-unitary evolution 
to be compatible with both energy-momentum conservation and locality. Moreover, 
in Sec. 5 we argue that energy-momentum conservation is not merely allowed, but 
even mandatory, for a large class of non-unitary systems with gravity. Finally, we 
apply these general results to the case of Hawking radiation in Sec. 6, and make some 
more general remarks in Sec. 7. 


2 A critical re-examination of BSP arguments 

In this section we re-examine some of the crucial arguments used in [13], with emphasis 
on those arguments which we find problematic. 

BSP find that a non-unitary evolution of the density matrix p can be described 
by a differential equation of the form 

p = N[ifo,d+ pQ f Q“ - 2 VVQ't, (1) 

which is Eq. (9) in their paper [13]. Here Q° are unspecified hermitian operators, h a p 
are unspecified c-numbers, and Hq is an operator which can explicitly be expressed 
in terms of Q a . (Actually, they make a typo by omitting the h a p factor in their 
Eq. (9), but it does not influence their further analysis.) Then, at the beginning of 
their Section 4, they make the following argument: 

“The failure of energy conservation can be seen from the following obser¬ 
vation: What if the theory did possess some hermitian operator H (not 
necessarily equal to H 0 ) which was conserved by the dynamics? Then any 
p which was a function only of H could not change under the action of 
(9). However, this is possible only if (9) contains only operators Q which 
are simultaneously diagonalizable with H. Unless H has highly degener¬ 
ate eigenvalues (a property which would exclude it as a good candidate for 
the energy), this is a serious restriction on h Q p, especially if Q a must be a 
local operator rather than a global charge.” (our italics) 
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After that argument, in the rest of their paper they study in detail these restrictions 
on h Q p and explain how these restrictions violate either locality or energy-momentum 
conservation. 

The italicized part of the quote above represents the core of our criticism of the 
BSP argument. Indeed, the italicized part shows that they are aware of a possible 
loophole for their argument; the argument is invalid if H has highly degenerate eigen¬ 
values. Nevertheless, the quote shows that they dismiss this loophole by saying that 
it “would exclude it as a good candidate for the energy”. Unfortunately, they do not 
explain why a Hamiltonian with highly degenerate eigenvalues would not be a good 
candidate for the energy and we cannot find any reasonable justification for such a 
claim. Just the opposite, we claim that a good candidate for the Hamiltonian with 
many degrees of freedom must have highly degenerate eigenvalues. 

The physical reasons for high degeneracy of the Hamiltonian are explained in 
detail in Sec. 3, but let us here present a simple mathematical argument. In a system 
with N degrees of freedom, the Hilbert space is usually spanned by a set of basis 
vectors of the form 

| 6 , 6 > ( 2 ) 

parameterized by N real (continuous or discrete) parameters £i, £2, • • ■, £jv- One of 
the parameters, say £ 1 , usually represents the energy E, so that all other parameters 
£2, ••■,£/v describe degeneracy for that energy. Hence, for large N, one expects a 
rather large degeneracy. 

Anyway, as is made clear in the quote above, the restrictions on h af 3 studied in [13] 
are only relevant when H is not highly degenerate. In other words, the whole BSP 
argument for violation of either locality or energy-momentum conservation rests on 
the assumption that the Hamiltonian is not highly degenerate. But this, we claim, is 
a false assumption, so their argument is flawed. Indeed, when H is highly degenerate, 
then there is a large set of mutually independent operators which commute with 
H. In this case, the analysis in [13] does not exclude the possibility that all the 
operators Q° in (1) belong to this large set, in which case the evolution (1) conserves 
energy without any restrictions on h Q p. We discuss such a conservation of energy and 
momentum in more detail in Sec. 4, by studying a more convenient Lindblad form of 
the evolution equation (1). 


3 Physical Hamiltonian is highly degenerate 

Why does the Hamiltonian must be degenerate? In essence, because no changes 
would ever happen in a universe with a conserved but non-degenerate Hamiltonian. 
If there was only one state with given initial energy, then energy conservation would 
forbid any transition to a new state different from the initial one. Nothing would ever 
happen in such a universe, in a clear contradiction with observations. 

But how much degenerate a Hamiltonian can be? Can the number of different 
states with the same energy be relatively small? In classical mechanics with N degrees 
of freedom, the total phase space is 2N- dimensional, while the phase subspace of a 
fixed energy E has dimension 2N—1. (Indeed, one of the authors of the BSP argument 
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considers this fact to be a part of the theoretical minimum one needs to know to start 
doing physics [21].) For a realistic system with large N, a (21V — l)-dimensional space 
of states with the same energy corresponds to a rather large classical degeneracy of 
the Hamiltonian. 

How about quantum systems? When the evolution is unitary, a convenient way to 
describe changes in quantum physics is to study the S'-matrix with matrix elements 
(f\S\i). As explained in any textbook on quantum held theory (including the ones 
[22, 23] by the other two authors of the BSP argument), the S-matrix conserves 
energy, i.e. states |i) and |/) are degenerate eigenstates of the Hamiltonian. Indeed, 
the number of states with the same energy is very large. For instance, states such as 

|w), 

all have the same energy E = co. (Here |ca) is a one-particle state with energy co, the 
state |^, f) is a two-particle state each of which has energy etc.) In a calculation 
of a total cross section, the integration over the quantum phase space of all states 
allowed by conservation laws may sometimes even lead to an IR divergence. 

Of course, the common energy E = co of the states in (3) is the energy of the 
free Hamiltonian. Can interactions reduce the degeneracy? In general, if |1) and |2) 
are two states with the same free energy Eq, the interaction may cause a split in the 
energy spectrum, so that |1) and |2) get different energies 

Ei = Eq + 5 Ei, E 2 = Eq + 5 E 2 , (4) 

respectively. Thus, if the measured system is in the superposition |1) + |2) before 
the interaction, during the interaction the measured system may end up in one of 
the states |1) or |2). However, such a transition must conserve energy of the whole 
system, which includes both the measured system and the source of interaction. If, 
for instance, during the interaction the measured system ends up in the state |1), then 
its change of energy by 5E\ must be compensated by a change —5E\ in the energy 
of the source, so that the total energy does not change. An analogous statement, 
of course, is also true when the measured system ends up in the state |2). Hence, 
even though the subsystem states |1) and |2) may have different energies due to 
interaction, the corresponding two states for the whole system must have the same 
energy. This demonstrates that the interaction does not reduce degeneracy of the 
total Hamiltonian. 
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4 Energy-momentum conservation in open quan¬ 
tum systems 

In general, energy and momentum of an open system do not need to be conserved; a 
part of energy and momentum can be transferred to the environment. Nevertheless, 
it does not mean that energy and momentum can not be conserved in an open system. 
In this section we explain in detail how a large degeneracy of the Hamiltonian implies 
a possibility of energy-momentum conservation, without violation of locality. 
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For simplicity we work in the Markov approximation, in which case, as derived 
in Appendix A, the density matrix p(t) of the open system satisfies the Lindblad 
equation [24] 

r l l 

P = -i[ H , P\ + L kP L l ~ ~ L k L kP - ^pL\L k . (5) 

k +0 L Z Z \ 

The Lindblad equation is known to be the most general Markovian and time-homogeneous 
differential equation for p(t) [24], Indeed, the Lindblad form (5) can be obtained 
from (1) by diagonalization of h a p [17]. The Lindblad equation is discussed in many 
monographs on quantum mechanics [25, 26, 27, 28, 29, 30] and even in some general 
quantum-mechanics textbooks [31, 32, 33]. 

The operators H, L k and p appearing in (5) act in a Hilbert space H. This means 
that the objects H , L k and p live in the product space H 2 = H® H*, where H* is 
the space dual to H. The Hilbert space H can be decomposed as 

H = (6) 

E 

where He is the Hilbert state of all eigenstates of the Hamiltonian H with the eigen¬ 
value E. In other words, all states in He have the same energy E , and the dimension 
of He measures the degeneracy of the Hamiltonian in the sector of energy E. The 
Hamiltonian can then be written in the spectral form 

H = Y,Y. e \E, i){E,Q, (7) 

E S 

where \E, £) are eigenstates of H, while £ is a collective label (compare with (2)) 
which labels different states with the same energy E. 

Now assume that p(0) G H 2 Eq , i.e. that initially the system has some definite 
energy E 0 . Then (7) implies [p(0),H] G H 2 Eq , implying that the first term in (5) 
conserves energy. Possible violation of energy-conservation can only arise from the 
second term in (5). However, if each L k lives in some subspace H%, or in other words, 
if for each L k there is energy E such that 

L k e He, ( 8 ) 

then also L\ G H 2 e and we have 

L kP ml e Keo for E = E 0 , 

L k p(0)L{ = 0 for E ± E 0 , (9) 

and similarly for other sub-terms in the second term in (5). This shows that the energy 
is conserved by (5) when L k satisfy (8). But the constraint (8) allows nontrivial L k 
only when the dimension of He is larger than one, i.e. only when the A-sector of the 
Hamiltonian is degenerate. Hence, when the Hamiltonian is degenerate, then L k can 
be such that (5) describes a non-unitary evolution which conserves energy. 

What about conservation of the 3-momentum P? We assume that the momentum 
is conserved by the purely Hamiltonian evolution, i.e. that 

[H, P] = 0. (10) 
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But then P has the same eigenstates as H and can be written in the spectral form 
similar to (7) 

p = EEp( s .fll £ .0< £ .«l. (n) 

E ( 

where p(E, £) are eigenstates of the operator P, i.e. P|if,£) = p(E,£)lE,£). Thus, 
denoting by H p the Hilbert space of all states with the same total momentum p, we 
see that L k can be such that they satisfy both (8) and 

L k e U\. ( 12 ) 

In this case (5) conserves both energy and momentum. Hence we see that a large de¬ 
generacy of the Hamiltonian may lead to conservation of both energy and momentum. 

In [13], where it was assumed that degeneracy is not present, it was shown that 
conservation of energy-momentum is incompatible with locality. Let us see how the 
presence of degeneracy avoids that conclusion, i.e. how both locality and energy- 
momentum conservation can be present in the presence of degeneracy. For that 
purpose, let us assume that the first sub-term in the second term in (5) has a local 
form 

J2 L kP L l = [ d 3 x£(x)p£ t (x), (13) 

k J 

and similarly for the other sub-terms. By introducing the Fourier transforms 

c ^=Swf e ~ l " L(p) ' £,(x)= /S eip,XiV)> (14) 

inserting them in (13), and using the identity 

J d 3 x e -hP-p')x = (2 tt) 3 5 3 ( P - p'), (15) 

one obtains 

J2 L kP L l = t 77^L(p)pL\p). (16) 

How do operators L(p) and Lqp) act on pi In the absence of degeneracy, one might 
argue that they change the total momentum of p by p, which violates the conservation 
of momentum. Indeed, this is essentially what BSP argued in [13] after their Eq. (21) 
(with an unimportant difference that they worked with operators <5"(p) instead of our 
Lindblad operators L(p)). However, in the presence of degeneracy such an argument 
is unjustified. When P is degenerated then L(p) can satisfy (12), in which case the 
momentum is conserved. Indeed, such operators L(p) can easily be constructed from 
any two different states |p, £) and |p, £') with the same momentum p, by taking 

L{ p) oc |p,f)(p,f'|- (17) 

Likewise, when H is degenerated then L( p) can satisfy (8), in which case the energy 
is conserved. In other words, the degeneracy implies that locality is compatible with 
energy-momentum conservation. 
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At this point it is also instructive to mention a related error made by Srednicki 
in [15]. He claimed that in quantum field theory there are very few operators avail¬ 
able which commute with the Hamiltonian, and that all of them are global (that is, 
integrals over all space of a local density). As the very few examples he had in mind, 
he mentioned Hamiltonian 77, momentum P, and a conserved charge. But contrary 
to his claim, we see that there are many different operators of the form (17) which 
commute with the Hamiltonian (7). They can be used to construct various local op¬ 
erators of the form (14) which also commute with the Hamiltonian. Finally, they can 
be integrated over space to construct various global operators which also commute 
with the Hamiltonian. The high degeneracy of 77 and P opens many possibilities 
for construction of various local or global operators which commute with 77 and P. 
Of course, unlike the few operators such as 77, P and the conserved charge, most 
of these operators cannot be easily expressed in terms of the canonical field opera¬ 
tors and their conjugate momenta. Nevertheless, all these operators share the same 
physical Hilbert space 77 on which they act. 


5 Open systems with gravity 

So far we have shown that energy-momentum in an open system can be conserved, 
but we have not claimed that it must be conserved. But when gravity is also taken 
into account, we can strengthen our claims. For a certain wide class of open systems 
involving gravity, in this section we argue that energy-momentum of the open system 
must necessarily be conserved. 

Let us start with energy-momentum conservation in classical general relativity. 
Let the metric g^{x) approach the Minkowski metric r) jW at infinity. By writing 
= r] fll/ + the Einstein equation can be written in an exact but non¬ 

covariant form [34] 


~ = —8vr G n t^(x). (18) 

Here ris a pseudo-tensor (i.e. a tensor under Lorentz transformations, but not 
under general coordinate transformations) which satisfies 

d u T^(x) = 0. (19) 

All the indices are raised and lowered by r) lxv and rR$ is the part of the Ricci 
tensor linear in h lw , and 


Tr(x)=Tp*°'(x)+T*r(x), ( 20 ) 

where 7)“ atter is the matter energy-momentum tensor and T^ av is the gravitational 
energy-momentum pseudo-tensor [34], Introducing the global pseudo-vectors 


P!‘ 


matter 


= [ d 3 xT^l tter (x), 


pR — 

1 grav 


<PxT^(x) 


( 21 ) 
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from (19) we see that the total energy-momentum pseudo-vector 

= Matter + P& av ( 22 ) 

is conserved. In addition, the Einstein equation (18) shows that the total energy- 
momentum can be interpreted as the source of gravity. 

Now consider a compact simply connected 3-dimensional region 1Z in space. The 
whole space consists of the interior of 1Z and the exterior of 1Z. Hence the total 
energy-momentum can be decomposed as 

P^PZ + PZt, (23) 

where P(f and P^ xt are the internal and external energy-momentum, defined by the 
corresponding internal and external regions of integration in (21). 

Now let the open system be defined as the collection of all degrees of freedom in the 
exterior of 1Z. The exterior can exchange energy-momentum with the interior, so P e xt 
does not need to be conserved. At first, one might interpret this as non-conservation 
of energy-momentum in the open system. However, the interior energy-momentum 
Ph is not really hidden to the exterior observer. The interior energy-momentum is 
the source of gravitational held, and this gravitational held produced by the interior 
energy-momentum can be measured even at the exterior. Hence, by measuring the 
external gravitational held, the external observer can determine Pjf. (This is similar to 
classical electrodynamics, where one can determine the interior charge by measuring 
the exterior electromagnetic held.) In this way the exterior observer can know not 
only P e xt, but also Pf. Consequently, the total energy-momentum (23) is known 
to the external observer, so effectively such an open system can be said to conserve 
the total energy-momentum. More precisely, the open system conserves information 
about the total energy-momentum, but operationally this can be thought of as being 
the same as conserving the total energy-momentum itself. 

Such an operational way of thinking in terms of information is particularly natural 
in the context of quantum- information [29]. Indeed, the density matrix p of the open 
quantum system describes the information available in the system. Hence, it seems 
reasonable to expect that in quantum gravity the external density matrix p conserves 
energy-momentum, at least at distances much larger than the Planck length, where 
the Einstein equation (18) is expected to be a good approximation. 

Let us express this expectation in a more quantitative form. Assume that initially 
the external state is a pure state |'0(O)) ex t with a definite total energy E. At a later 
time t the external system can become entangled with the internal one, so the total 
state becomes a superposition of the form 

l*«> = E£CrfWMiJs,0™., (24) 

where C x ^(t) are some suitably normalized coefficients. Due to the conservation of 
total energy in the exterior, all the external states | E, £) ext have the same total energy 
E equal to the initial exterior energy. (For simplicity, the conservation of 3-momentum 



is not considered explicitly, but it is not difficult to include it as well.) These external 
states can be further decomposed into the gravitational and matter parts as 


l-^oOext | ^grav) Cgrav) | ^matter ? ^matter)) (25) 

where E gl . av = E — E mattei , so (24) can be written in a more detailed form 

l«V)> = EEEE "X, -^matter£grav^matter (^0 I in 

X ^matter £grav ^matter 

®| E- -^matter) Cgrav) [-^matter) 'Cmatter) • (26) 


The external state pit) is obtained by tracing the pure state |T(f))(T(f)| over the 
internal degrees of freedom, so (24) leads to 

pit) = J2 w dt) \ E 5 £ ) ext ext <£.a (27) 

s 

where w^(t) = \ E xt(t)\ 2 - The state (27) is a mixed state with a definite energy E. 

Recalling that the initial external state was the pure state p(0) = |V ; (0)) e xt ext ("0(0)1 
with the same energy E, we see that the evolution (27) is a non-unitary evolution 
which conserves energy. 


6 Application to Hawking radiation 

Now let us apply the general results of the preceding sections to the case of Hawking 
radiation [1], The interior region 7 Z is the interior of a black hole, which is a special 
case of the situation studied in Sec. 5. For simplicity, we consider only spherically 
symmetric modes of radiation, so the conservation of 3-momentum is trivial. Hence 
only the energy conservation is considered explicitly. 

6.1 Energy conservation by Hawking radiation 

Let the initial exterior state be a pure state 

WO)) = |M)|0). (28) 

Here | M) is the external state of gravity the source of which is a black hole with 
initial mass M, while |0) is the matter vacuum. After some time At, the total state 
takes the form [35] 

I*(a*)> = EE ^^radCrad (At) | - -®rad? £rad)in|Af E iad) \ E ra,di £rad); (29) 

-^rad £rad 

which is a special case of (26). Here | M — 77 rad ) is the external gravitational state cor¬ 
responding to the black-hole mass M' = M — E r&( \ , while I7 rad is the energy of Hawking 
radiation in the state | A rad , £™d) • The gravitational state | M — A rad ) does not have 
any degeneracy parameter £ grav , which is a consequence of the no-hair theorem. (This 
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theorem [36, 37] states that the gravitational held external to the stationary black 
hole reveals no other information about the black hole except its mass, charge and 
angular momentum, and for simplicity we assume that charge and angular momen¬ 
tum are zero.) The external states of Hawking radiation \E rad) £ ra d) and the internal 
states | — .£/ ra d, £ ra d)m are parameterized by the same parameters E ra d and £ ra d, which 
corresponds to the perfect correlation of external Hawking radiation with internal 
Hawking quanta. (Of course, the internal Hawking quanta are not the only quantum 
degrees in the interior. However, the other internal degrees will not play explicit role 
in our calculations, so we suppress them.) The Hawking-radiation states \E rad , 6-ad) 
with the same E ra d but different £ ra d include different states of the form (3). 

The Hawking radiation has a thermal spectrum, which means that 


g—/t^rad 

|Ce.^^(A «)| 2 = z(At)' 


(30) 


Here Z(At) is a normalization constant to be discussed later and f3 1 
Hawking temperature [1] 


T 


1 

87rM’ 


T is the 
(31) 


in the natural units Ti = c = Gn = k-Q = 1. The external state associated with (29) is 
the mixed thermal state 


p(At) = E E l^ad^d(A*)| 2 1 M,E rad j £rad) (M, T/ ra d, £rad 15 (32) 

^rad £rad 

where 

| M, E Iad , £ r ad) = \M — E rSLd )\E iad , ^ rad ) (33) 

are states all having the same total external energy M. 

6.2 Violation of unitarity 

The evolution (29) is unitary, which is related to the fact that it is obtained with the 
aid of a Bogoliubov transformation [38]. However, the creation of internal Hawking 
quanta | — E ra d,£ ra d)m increases the internal entanglement entropy of the black hole. 
In this way the black hole entropy can easily become larger than the Bekenstein- 
Hawking entropy S = Aj 4, where A = 4rcR 2 is the area of the black-hole horizon 
at the radius R = 2 M. This is a problem, not only because there are many reasons 
to believe that S = A /4 should be the correct black-hole entropy, but even worse, 
because in this way the black-hole evaporation could eventually create a very small 
and light black hole with an arbitrarily large amount of entropy. This is one of the 
ways to formulate the infamous black-hole information paradox. 

This form of the paradox can be removed by assuming that the internal Hawking 
quanta are destroyed by a non-unitary process near the black-hole singularity. But 
then the external mixed state (32) becomes a state of a closed system. That means 
that a closed system evolves non-unitarily from a pure to a mixed state, which the 
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standard quantum theory does not allow. This represents another form of the black- 
hole information paradox. 

Of course, it is conceivable that unitarity of standard quantum theory might 
break down at black-hole singularities. But, as said in the Introduction, one of the 
most frequent arguments against non-unitary evolution for the closed system is the 
BSP argument [13] that such evolution would be incompatible with either energy- 
momentum conservation or locality. Yet, in this paper we have shown that such an 
argument is incorrect. Indeed, the non-unitary evolution (32), now interpreted as 
the evolution for a closed system, conserves energy-momentum. On the other hand, 
the creation of Hawking particles at the horizon is a local process [39, 40, 41]. If, 
as we suggested above, the true violation of unitarity happens near the black-hole 
singularity where the internal Hawking quanta are destroyed, then the violation of 
unitarity does not seem to violate locality, or at least not at distances much larger 
than the Planck length. 

6.3 Lindblad operators for Hawking radiation 

The general argument in Sec. 4 for compatibility between non-unitarity and energy- 
momentum conservation involved the Lindblad equation. To see the connection be¬ 
tween (32) and the Lindblad equation, let us now calculate the Lindblad operators 


corresponding to (32). 

For our case, Eq. (49) in Appendix A can be written as 

in( — A rad , £ rad |\I/(At)} = A^ rad £ rad (At) | M) |0) . (34) 

Using (29), this gives 

AE rad £ rad (At)| M) |0) = C Eiad z ia d (At) \M - A rad ) |U rad , £ rad ), (35) 

which in the notation of (33) can be written as 

Ae rad £ rad (At)|M, 0,0) = U Erad £ rad (At)|M, A rad , £ rad ). (36) 

Hence the Kraus operators A Ei ^ rad (At) are 

^E rad $ rad (At) = C Eiad ^ iad (At)\M, £ rad , £ rad ){M, 0,0|. (37) 


In (29), |C(E rad £ rad (At)| 2 is the probability that the system will be found in the 
corresponding state after the time At. For a relatively short time At (but not shorter 
than the Zeno time [42]), this probability increases linearly with time. Therefore we 
can write 

|As rad q ad (At)| 2 = At|C^ rad £ rad | 2 ,. (38) 

where |A^ rad ^ ad | 2 is the transition probability per unit time which does not depend 
on At. Hence (30) and (38) imply 

p /^-^rad 

l^kadh-aJ 2 = Y' ’ ^ 
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where Z' — At Z(At) does not depend on At. 

To find the constant Z r , we calculate the total average power of Hawking radiation 

= £ E^lc^J 2 = Y' £ £«-«-*-■ («) 

-^rad £rad ^rad £rad 



The left-hand side of (40) can also be calculated from the Stefan-Boltzmann law of 
black-body radiation 

(^)=«A T \ (41) 

where o = 7r 2 /60 is the Stefan-Boltzmann constant, A = AnR 2 = 16nAL 2 is the 
black-hole area, and T is the Hawking temperature (31). This gives 


/ dE ISU d \ CT 

\ dt / - M 2 ’ 


(42) 


where o' = a/2567r 3 . Hence (40) implies that the normalization constant Z' is 

v/ _ E Srad E ?rad ^rade-^ 

z - VJm 2 ' ( } 

Now the Lindblad operators can be calculated by applying the general equation 
(52) in Appendix A to our specific equation (37). Using also (38) and (39), this finally 
gives 

e -/3E lad /2 

^radSrad = ^ ^rad, frad) ( M, 0, 0|. (44) 

Note that the phases of the Lindblad operators are irrelevant, because the phases 
cancel out in (5), in products involving both Lk and L\. As the states \M, U ra( j, ^ ra d) 
and \M, 0,0) both have the same total energy M (recall Eq. (33)), we see that the 
Lindblad operators (44) conserve energy. (The conservation of 3-momentum is trivial, 
as explained in the introductory paragraph of Sec. 6.) 

We also see that the Lindblad operators (44) depend on the black-hole mass M. 
Therefore, they can be treated as time-independent operators as long as the black- 
hole mass does not change much. This is a very good approximation for quite a long 
time, because (42) implies that mass of the black hole changes significantly only after 
a time of the order of M 3 . Indeed, denoting by M(t ) the average black-hole mass at 
time t, we see that (42) and the conservation of total average energy imply 

dM (t) o' 

dt Jl 2 (t) 

The solution of this differential equation is 

M(t) = ^M 3 - 3 o>t, 

where now the initial black-hole mass is denoted by Mq. 


(45) 

(46) 


12 



In this way the Lindblad equation (5) can be improved by replacing the time- 
independent Lindblad operators L*. with the time-dependent ones L^it). The time- 
dependent Lindblad operators are given by (44) in which M is replaced by M(t). 
Of course, this is still an approximation, but now a better one. In particular, the 
resulting master equation is still Markovian, but not longer time-homogeneous. In 
addition, the resulting time-dependent Lindblad operators still conserve energy and 
momentum. 


7 Discussion and outlook 

In this paper we have argued that one of the main alleged problems with non-unitary 
Hawking radiation, the problem of violation of either locality or conservation of 
energy-momentum, does not really exist. Owing to a high degeneracy of the Hamil¬ 
tonian and momentum operators, non-unitary Hawking radiation is fully compatible 
with both locality and energy-momentum conservation. This eliminates one impor¬ 
tant obstacle for the possibility that Hawking radiation might really violate unitarity. 

This, however, does not imply that unitarity violation by Hawking radiation is 
without problems. There are also other good reasons (not eliminated by results of 
the present paper) to think that evolution should be unitary. 

The first reason is the fact that known non-gravitational quantum theories are 
unitary, suggesting that quantum gravity should not be different. Nevertheless, grav¬ 
ity is different by being a theory of space and time. In particular, unitary evolution 
requires a globally hyperbolic spacetime, while general relativity contains solutions 
which are not globally hyperbolic. This suggests the possibility that, in the presence 
of gravity, quantum theory should be generalized to a theory in which the notion of 
unitarity does not need to refer to a unitary time evolution [43, 44, 45, 46]. In this 
way the non-unitary time evolution of Hawking radiation may not be in contradiction 
with a generalized concept of unitarity. 

The second reason, especially popular in the string-theory community, is the 
AdS/CFT correspondence [47], or more generally the gauge/gravity correspondence 
[48]. According to this correspondence, the quantum theory of gravity is mathe¬ 
matically equivalent to a quantum held theory on the boundary. Since the latter is 
known to be unitary, it follows that the former must also be. This, indeed, is a very 
strong argument for unitarity of Hawking radiation, but is not without possible loop¬ 
holes. First, many aspects of this correspondence depend on the assumption of string 
theory, while we still do not have an experimental proof that string theory is the 
correct theory of Nature. Second, many aspects of the correspondence depend on the 
assumption of supersymmetry and AdS background, while the Universe we observe 
does not look supersymmetric and we do not seem to live on an AdS background. 
In the absence of supersymmetry and AdS background, the present evidence for the 
correspondence is still incomplete. Third, even if these problems are merely technical 
and the correspondence is exact even without supersymmetry and AdS background, 
there is an intriguing suggestion that the field theory part of the correspondence may 
in fact be a non -unitary gauge theory [49]. 
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To conclude, the question of unitarity of Hawking radiation seems to be an open 
question. In our opinion, non-unitary Hawking radiation is still a viable possibility 
worthwhile of further investigation. 
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A Lindblad equation 

Since many intended readers of this paper are probably not familiar with the Lindblad 
equation, for the sake of completeness we present here some basics, by giving simple 
derivations and explanations of the relevant equations used in the rest of the paper. 
Our derivations are to a large extent based on [33], were some further illuminating 
details can be found. 

Let l'L(t)) be the full time-dependent state of a closed system, and let the system 
be split into two subsystems, one of which is observed and the other unobserved. Let 
| k) be a basis for the unobserved subsystem. Then the density matrix of the observed 
subsystem is 

P(t) = TV mob ,|>I>(«))((P(i)| = £<fc|¥(«)>(t(S)|A : >. (47) 

k 

This can be written as 

p{t) = x^k(*)p(°)4to> ( 4s ) 

k 

where A^it) are the so-called Kraus operators which act only in the Hilbert space of 
the observed subsystem. In particular, if the observed subsystem is initially in the 
pure state p(0) = 1-0(0)) (0(0)1, from (47) and (48) one can see that the action of 
Kraus operators is given by 


«*(«)> = /l*(«M0)>. (49) 

In general, the density matrix must satisfy Trp(t) = 1 for any p(0), so (48) leads to 
the constraint 

E4(l)4(i) = l. (50) 

k 

Now consider p(At) at a small time At, so that 

p(At) = p(0) + A p. (51) 

Assuming that one of the Kraus operators (call it A 0 ) is close to the unit operator, 
we can write 


Ag(At) — 1 + At (Lq — iH), 

kb,(At) = \ZAtLfc for k ^ 0, (52) 
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where L 0 and H are some hermitian operators. In this way, in the lowest order in 
At, (48) leads to the equation 


p{ 0) + A p = p( 0) + At {L 0 , p(0)} - i[H, p(0)] 


E Lk P m\ 


(53) 




Strictly speaking, this equation is correct only initially, i.e. for p(0). Nevertheless, 
if the density matrix p{t) quickly “forgets” its history, or equivalently, if the time 
evolution of the observed subsystem is sufficiently slow so that there is enough time 
to “forget” the history, then one can use the Markov approximation, i.e. one can 
consider (53) to be approximately valid for p(t ) at any t. 

Or let us be slightly more precise [33]. When the unobserved degrees interact 
with the observed ones, the interaction usually creates entanglement between them, 
which makes p mixed. However if, after the creation of entanglement, the unobserved 
degrees are removed from the vicinity of the observed ones, then the unobserved 
degrees have no further influence on the evolution of the observed degrees. In this 
way the density matrix p at a given time t effectively “forgets” that it was mixed 
due to entanglement with the unobserved degrees, and only “remembers” that it is 
now mixed. In other words, all the information relevant for further evolution of p is 
contained in the p{t) itself, not in its environment. 

In this approximation, and in the limit At —> 0, (53) becomes the differential 
equation 

= {L 0 , p(t)} - i[H, p(t )] + L kp{t)L\- (54) 

ai o 

Furthermore, the condition Tr p(t) = 1 implies 


Tr A!! = o, 

dt 

so from (54) one finds a constraint similar to (50) 

2Fo + Y, L\L k = 1. 

k^ 0 


(55) 


(56) 


This can be used to eliminate Lq from (54), so (54) can be written in the final form 
dp 


dt 


= —i [H, p\ + Y 


M o 


LkpL\ — -L\L k p — -pL\L k 


(57) 


This is a master equation known as the Lindblad equation [24], while the operators 
L k are called Lindblad operators. 

We have not yet explained the physical interpretation of the operator H. All we 
know at this moment is that it is a hermitian operator acting in the Hilbert space 
of the observed subsystem. In general, the evolution (57) is not unitary. But in 
the special case in which all L k = 0 for k ^ 0, (57) reduces to the familiar unitary 
evolution generated by H. This shows that H can be interpreted as the Hamiltonian 
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for the observed subsystem. However, H depends also on the interaction with the 
unobserved subsystem, meaning that H is a modified Hamiltonian that includes the 
unitary influences of the unobserved subsystem to the observed one. For example, 
when a laser beam is reflected at a heavy mirror, it is a unitary process for the beam 
[33]. So if p is the state of the laser beam, then H is an operator acting only in the 
Hilbert space of the beam, in which case H contains both the free Hamiltonian of the 
beam and a term describing the influence of the mirror. 
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